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ABSTRACT 

It is generally accepted that magnetic fields generated in the nonlinear development of the transverse Weibel 
instability provide effective collisionality in unmagnetized collisionless shocks. Recently, extensive two and 
three dimensional simulations improved our understanding of the growth and saturation of the instability in 
colliding plasma shells. However, the steady-state structure of the shock wave transition layers remains poorly 
understood. We use basic physical considerations and order-of-magnitude arguments to study the steady state 
structure in relativistic unmagnetized collisionless shocks in pair plasmas. The shock contains an electrostatic 
layer resulting from the formation of stationary, magnetically-focused current filaments. The filaments form 
where the cold upstream plasma and the counterstreaming thermal plasma interpenetrate. The filaments are not 
entirely neutral and strong electrostatic fields are present. Most of the downstream particles cannot cross this 
layer into the upstream because they are trapped by the electrostatic field. We identify the critical location in 
the shock transition layer where the electromagnetic field ceases to be static. At this location, the degree of 
charge separation in the filaments reaches a maximum value, the current inside the filaments comes close to 
the Alfven limit, and the phase space distribution function starts to isotropize. We argue that the radius of the 
current filaments upstream of the critical location is about twice the upstream plasma skin depth. Finally, we 
show that some downstream particles cross the electrostatic layer and run ahead of the shock into the preshock 
medium without causing instability. These particles may play an important role in particle acceleration. 

Subject headings: instabilities — magnetic fields — plasmas — shock waves 



1. INTRODUCTION 

Relativistic collisionless shocks should accompany all rel- 
ativistic outflows, such as those associated with active galac- 
tic nuclei (AGN), microquasars, gamma-ray bursts (GRBs), 
giant flares from soft gamma-ray repeaters, and pulsar wind 
nebulae. Nonrelativistic collisionless shocks are also believed 
to be present in cosmological structure formation, supernova 
remnants, and planetary bow shocks. The composition of the 
relativistic outflow (e.g., e~e'^, e~p, or Poynting flux) is not 
securely determined in any of the relativistic sources, and the 
level of magnetization is poorly constrained. In some systems 
(e.g., AGN) it seems that the magnetic fields in the upstream 
medium of the shock are dynamically important, while in oth- 
ers (e.g., "external" GRB shocks), the upstream is effectively 
unmagnetized. In the latter case, the magnetic fields that are 
invoked to explain the observed synchrotron emission seem to 
be generated in the shock itself. 

Magnetic fields with energy density near equipartition with 
the kinetic energy density of the plasma can be generated by 
the nonlinear development of the transverse Weibel instabil- 
ity. The Weibel instability (.Weibel. 1959) arises when the ve- 
locity distribution of the plasma particles is anisotropic (e.g., 
lYoon & Davidson"198'A. The transverse Weibel instability 
(e.g., iFried 1959; Silva et al. 2002; Bret. Firpo. & Deutscl 
|2()04, 2005) results in the breakdown of the plasma into cur- 
rent filaments ( Lee & Lamp e 1973 ). Strong anisotropy arises 
in collisions in which charged particles interpenetrate. This 
has naturally led to suggestions that the transverse Weibel in- 
stability is responsible for the gro wth of magnetic fields in 
unmagnetized collisionless shocks ( Medvedev & Loeb 
IGru zinov 2001J. In this picture, magnetic fields generated 
in the shocks provide effective colUsionality in the otherwise 



collisionless transitional layer of the shock. The magnetic 
fields may also be responsible for diffusive particle acceler- 
ation in the shock. 

The growth and saturation of the transverse Weibel insta- 
bility in colliding cold unmagnetized plasma shells (both 
e'e"*" and e~p) has been extensively studied using three- 
dimensi onal particle-in - cell (P I Q simulations ( Silva et aD 
(2003li 'Fonseca et al." '2003t iFrederik sen et al. 20041 
Jaroschek. Lesch. & Treumann 120041 12005: Nishikaw a et alJ 
2005; Spitkovsky & Arons 2005, in preparation). These sim- 
ulations are computationally demanding; they must resolve 
multiple plasma skin depths in the direction perpendicular to 
shock propagation (henceforth, perpendicular direction), and 
a very large number of plasma skin depths in the direction 
parallel to shock propagation (henceforth, parallel direction) 
over which field growth and plasma isotropization and 
thermalization take place. PIC simulations have not yet 
converged to reveal the steady-state structure of the transition 
layer that persists after the transients associated with the 
initial collision of plasma shells have dissipated. 

Simulations of cold shell collisions show that c urrent fila- 
ments form parallel to the axis of collision (e.g., Silv a et all 
■200 3; Frede riksen et al. 2004) . The direction of current flow 
in neighboring filaments alternates. The currents are main- 
tained by the pinching force of the self-magnetic field. The 
maximum current that can flow through any given filament, 
precluding unphysical velocity-space configuration, is the 
Alfven limiting current (Alfven 1939; Davidson 1974; Hond^ 
I2OOO). The development of the transverse Weibel instability 
in cold, colliding shells involves the merging of current fil- 
aments ( Medvedev et al. 2005) and saturation of the Alfven 
Hmit (IKatnl200.5l) . 
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It has long been known that it is possible to construct ex- 
act, albeit unstable, solutions for collisionless shock transi- 
tion layers in one dimension, where kinetic energy of the up- 
stream plasma is converted into electrostatic potential energy 
(e.g., of plasma oscillations) as the plasma crosse s the transi- 
tion layer (Moiseev & Sagdeev 1963; Montgo mery & .Tovcd 
[l969). These solutions invoke particle reflection by the sta- 
tionary electrostatic potential of the transition layer and be- 
long to the class of "double layers." A double layer is an elec- 
trostatic structure which can appear within a current carrying 
plasma and sustain a significant net potential differenc e asso- 
ciated with at least two layers of opposite net charge iRaadul 
1 1989* and references therein). 

Moiseev & Sagdeev (1963) noticed that the interpenetra- 
tion of the reflected and the forward streaming plasma is un- 
stable to a Weibel-like mode. They conjectured that the scat- 
tering of particles by the electromagnetic fields generated in 
the instability gives rise to an effective viscosity that facili- 
tates momentum diffusion in the transition layer. Albeit pre- 
scient, this model has not developed into a self-consistent pic- 
ture of the shock transition layer 

Here we suggest a self-consistent steady state structure of 
a relativistic shock that propagates into unmagnetized pair 
plasma. In §|2j we sketch out a qualitative physical picture of 
the shock transition layer. In § |3] we introduce the reference 
frame and notation that is used in the rest of the paper. In §|3 
we discuss the motion of particles in the quasistatic electro- 
magnetic field of current filaments that form in the transition 
layer In § [S] we study the mechanics of the separation of the 
upstream flow into filaments of oppositely charged particles. 
In § |6l we study the kinematics of particles trapped by the 
electrostatic field in the transition layer, and show how these 
particles facilitate the onset of charge separation. In § we 
argue that the radius of the current filaments must be about 
equal to the plasma skin depth. In §|8] we show that particles 
in the high-energy tail of the thermalized downstream can es- 
cape into the upstream domain. In § |9j we discuss the com- 
pression and convergence to the hydrodynamic jump condi- 
tions, and we also speculate on the applicability of our model 
to e~p shocks. Finally, in §[^1 we summarize our main con- 
clusions. 

2. PHYSICAL PICTURE 

We consider an ultrarelativistic shock wave propagating 
into unmagnetized pair plasma. The structure of the shock 
transition layer is steady in the rest frame of the shock. In this 
frame the cold, upstream plasma arrives on one side of the 
transition layer, and thermalized downstream plasma leaves 
on the other side. Since the plasma is collisionless by assump- 
tion, electromagnetic fields inside the transition layer deflect, 
scatter, and redistribute the kinetic energies of the particles 
arriving from the upstream. Steady currents are present to 
support a steady electromagnetic field. 

The saturated state of the transverse Weibel instability, as 
seen in the simulations, consists of magnetically self -pinched 
current filaments. Initial filaments are produced at the point 
where the upstream encounters the transition layer These fil- 
aments merge to form structure on ever larger scales. We 
here focus on the first generation of filaments only, because 
it is key to understanding the subsequent evolution of the 
shocked plasma. The filaments form as the magnetic field 
displaces the cold upstream particles transversely, perpendic- 
ular to the direction of sh ock propagation (see, e.g.. Fig. 1 in 
iMedvedev &LoehlfT999l) . The filaments are in place when 



the displacement becomes similar to the separation between 
filaments. The filaments live quasi-statically for some dis- 
tance, and then become time dependent, as the electromag- 
netic fields turn from quasi-static to dynamic. 

The transverse Weibel instability requires interpenetrating, 
mutually counterstreaming plasma components. One of the 
components is the cold, upstream plasma. The other compo- 
nent is the downstream plasma, i.e., the particles originating 
in the upstream that are reflected in the downstream or in the 
transition layer and are now moving backward, toward the up- 
stream in the shock frame. Hereafter we refer to the latter 
particles as the "counterstream." The transition layer is where 
the upstream and the counterstream interpenetrate. 

The small-scale toroidal self-magnetic field of the current 
filaments cannot keep the counterstream particles trapped in 
the transition layer This is because the field configuration 
is open-ended facing the upstream and any particles affected 
by the magnetic field alone would escape by bouncing be- 
tween the magnetic walls until they have reached the open- 
ing. However, most of the counterstreaming particles have 
to be prevented from escaping into the upstream if hydrody- 
namic jump conditions are to be reproduced. We show that a 
parallel electrostatic field is present in the first generation of 
current filaments and that this field can contain the counter- 
stream. 

Since the filaments are quasi-static, the electric field can be 
derived from a potential associated with a quasi-static distri- 
bution of charge. The sign of the charge in neighboring fila- 
ments alternates, ensuring overall charge neutrality. Parallel 
electric field arises from parallel charge gradients within the 
filaments. The density of counterstreaming particles within 
the filament is less than that of the oppositely charged for- 
ward streaming particles, implying non-neutrality. Net charge 
does not prevent filamentation because only a small amount 
of counterstreaming particles is sufficient to ensure that mag- 
netic pinching can overcome electrostatic repulsion and sep- 
arate the oncoming upstream particles based on their charge. 
We show that the degree of charge-separation reaches maxi- 
mum at about the same point where the current reaches the 
Alfven limit. In a self-consistent picture, the radius of the 
first-generation filaments is about equal to the plasma skin 
depth of the upstream. 

We refer to the region upstream of the point of maximum 
charge separation as the "charge separation layer." In this pa- 
per we investigate the structure and properties of the charge 
separation layer in detail. The layer is located at the leading 
edge of the shock before most of the density compression. 
We present a picture in which the charge separation layer is 
described by an electrostatic double-layer that forms inside 
each magnetically-pinched current filament. The layer there- 
fore functions as a diode protecting the upstream from the 
downstream. The filaments themselves are a stationary end- 
product of the nonlinear transverse Weibel instability. 

Figure^schematically depicts the structure of a current fil- 
ament in the charge separation layer. The upstream particles 
are separated into filaments with a typical radius A of the or- 
der of the upstream plasma skin depth. We set the electrostatic 
potential to zero at x = -oo (far left in the figure). The poten- 
tial is also zero at the midpoint between filaments (the top and 
bottom dot-dashed line in the figure). A filament of positive 
current, carried mainly by the positrons from the upstream, 
produces a potential ridge (central dot-dashed line), while a 
filament of negative current produces a potential valley. As 
the magnetic field displaces an upstream positron toward the 
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Fig. 1. — Schematic representation of a positron within the charge sepa- 
ration layer; an electron filament would equivalent with the direction of the 
magnetic field B, electric field E. and current I, reversed. The conical en- 
velope {solid lines) is the trajectory of the last electron to leave the growing 
positively-charged domain {solid ellipse). The top and bottom dot-dashed 
lines separate the positron filament from the neighboring electron filaments. 
The central dot-dashed line is the filament center. We also show an equipo- 
tential surface (dashed lines). 



positive filament, the positron climbs the electrostatic ridge. 
An electron located in what is to become a positron filament 
is displaced out (solid line in the figure). 

The absolute height of the electrostatic potential at the point 
of maximum charge separation (xq in Fig.Q is comparable to 
the kinetic energy of an upstream particle. This is required 
for the counterstream to be decelerated and reversed, as the 
kinetic energy of the downstream particles, among which the 
counterstream particles originate, is smaller by a factor of 
~ 3\/2 than that of the upstream particles. The same poten- 
tial decelerates the upstream particles as they cross the charge 
separation layer (see Fig. |2j- If the upstream flow remains 
relativistic as it traverses the layer, compression must occur in 
the downstream, where time-dependent electromagnetic fields 
can scatter and isotropize the particles. 

To be contained within the transition layer of the shock, 
most of the counterstream must already be charge-separated 
as it enters the charge separation layer at xq. Namely, a 
counterstream positron must be contained within the poten- 
tial valley of the filament carrying upstream electrons, and 
vice versa (an equipotential surface is outlined with dashed 
lines in Fig.Q. As the positron attempts to climb out of the 
potential valley of the electron filament, it reverses direction 
and proceeds back toward the downstream (see Fig.|3}- 

If the energy distribution of the counterstream particles is 
thermal at xq, a small fraction of the particles in the exponen- 
tial tail have sufficient energy to escape the transition layer 
and run ahead of the shock. We show that plasma stability 
tightly limits the flux of counterstreaming particles that can 
run ahead of the transition layer into the untouched upstream. 
These particles may play important role in diffusive particle 
acceleration in shock waves. 

3. NOTATION 

We work in the reference frame in which the shock is at 
rest. Unless otherwise stated, all densities, velocities, mo- 
menta, and fields are written in this frame. It is common in the 
literature, however, to write particle densities in the local rest 
frame of the plasma ("fluid frame"), a nd to write momenta 
in the rest frame of the upstream (e.g., iBlandford & McKe3 



Table 1 . Notation 



Quantity 


Our Notation 


Ultrarelativistic Limit 


upstream in shock frame 




\/272 


downstream in s. f. 


7+CX) 


3^2/4 


upstream density in s. f. 


n-oo 


\/272>ii 


downstream density in s. f. 


n+oo 


3\/272(ii 


compression factor in s. f. 
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119761) . Lorentz factor of the unshocked upstream medium in 
the rest frame of the shock, which is usually denoted by F, we 
denote with 7-00- Similarly, we denote the Lorentz factor of 
the downstream with 7+00 • The densities of the upstream and 
downstream medium far from the transition layer are denoted 
with n_oo and n+00, respectively. 

In Table [U we express these quantities in terms of the no- 
tation of Blandford & McKee ( 1976), where 72 is the Lorentz 
factor of the downstream in the rest frame of the upstream; 
til and n2 are the rest frame densities of the upstream and 
the downstream, respectively, and 72 is the ratio of specific 
heats. Note that hydrodynamical jump conditions dictate that 
n+00 /n-oo = 3. All quantities referring to the counterstream 
particles are marked with a prime. Let xq denote the point of 
maximum charge separation; the quantities measured at xq are 
marked with subscript zero. 

4. BALLISTICS 

We here consider the ballistics of a particle moving inside a 
current filament. We use cylindrical coordinates (x, r, 9) where 
X is the symmetry axis of the filament parallel to the direc- 
tion of shock propagation. Since counterstream particles have 
large perpendicular velocities (Pr ~ 1), they cross the mid- 
dle of the filament many times. Therefore, the quantities per- 
taining to the counterstream particles can be evaluated at the 
center of the filament r = 0, e.g., /3'(x) = /3'(x,Q), where the 
point of maximum charge separation xq is where an average 
particle has executed one oscillation toward the center of a fil- 
ament. The upstream particles on average cross the filament 
only once between x = -00 and xq and thus their radial depen- 
dence must be retained. 

We approximate a filament with a cylinder of radius A. The 
filament touches neighboring filaments of the same radius that 
carry equal and opposite current. In reality, a filament inter- 
faces with the neighboring filaments over its entire area in 
such a way that no empty space remains between. Cylindrical 
filaments are thus an idealization. We assume that the prop- 
erties of the filament (net current, charge density, etc.) vary 
along the x axis on scales much larger than A, implying that 
the filament is locally two-dimensional. 

Because of the quasi-static nature of the filaments, the total 
energy of a particle, jmc^ + q(j>, is conserved. Here, q is the 
charge of the particle, 7 is the Lorentz factor, m is the mass 
of the particle, and is the electrostatic potential. Conserva- 
tion of energy for a particle originating in the upstream then 
implies 

7-oo'«c^ = j(t)mc^ + q4)(x(t), r(t)). (1) 

A cold upstream particle is deflected slowly towards the cen- 
ter of its filament. The equation of motion for the perpendic- 
ular momentum of the particle reads 

^ = q{E,- l3,Be), (2) 
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which, using p,- = ^mcdr/dt and d/dt = cfixd/dx, can be 
rewritten as 

mc^/3x^ [iPx^] = q(Er-Me). (3) 
dx \ dx J 

Note that for the current to pinch, the condition PxBg > E,- 
must be satisfied, implying that there exists a local reference 
frame in which the filament is charge-neutral and Er vanishes. 
We call this frame "the frame of the current." Let the velocity 
relating the frame of the shock to the frame of the current be 
P = -Er/Bg, and the corresponding Lorentz factor 7. 

A counterstream particle is electrostatically confined to the 
filament if (70- l)mc^ < -q'4)Q, where 70 is the Lorentz factor 
of the counterstream particle at (x, r) = (xo,0) and 00 is the po- 
tential at this location. A particle that is electrostatically con- 
fined to its filament is also bound to the transition layer and it 
must decelerate and reverse its direction (see Fig.|3}- A par- 
ticle which is not confined electrostatically may still be con- 
fined magnetically to the filament and reversed if /3^q ^ /3^q. 
A counterstream particle that is not confined to the filament 
runs away to x = -00. 

Near xq, charge separation is almost complete and the elec- 
tromagnetic field varies weakly with radius. We consider an 
upstream particle that is moving ultrarelativistically inside a 
filament in the same sense as the current in the filament, i.e., 
f3/3 > and 7^7^ 1 . The magnetic field tends to confine 
the particle within the current filament, while the electric field 
tends to repel it out of the filament. Given that \Er\ < \Bg\, 
the magnetic field wins by a small margin, and confinement 
is possible provided that the particle's momentum is not too 
large. It is straightforward to show that the particle remains 
confined provided that < XeBg/'-f'y^mc^. 

5. CHARGE SEPARATION 

We here derive an equation governing the degree to which 
the charges in the upstream are consigned to disjoint current 
filaments. Since the upstream plasma is neutral in the limit 
X -00, the region r < A in which the pinching of charges q 
takes place must be evacuated of charges -q. Here, A is the 
radius of the filament. Charge separation takes place gradu- 
ally. We trace the orbit (x(t),r{t)) of the last of the charges 
^ = to leave the region. This particle had moved just along 
the jc-axis before approaching the charge separation layer (see 
Fig. [3- The region r < r{t) can then be considered "charge- 
separated," while the region r < r{t) < A is still neutral. 

We make the simplifying assumption that the particle den- 
sity n and mean axial velocity of particles in the upstream 
is independent of r. The upstream plasma at x = -00 is undis- 
turbed, neutral, and cold, with density n_oo moving toward 
the shock with the Lorentz factor 7-00 in the frame of the 
shock. In this section we assume that the electromagnetic 
field is static, and thus that the electric field is irrotational and 
derived from a potential 0(x, r) which vanishes at x = -00. 

The electromagnetic field of neighboring filaments cancels. 
Therefore the dominant contribution to the electromagnetic 
force comes from a segment of length Ax ~ 2A of the closest 
(host) filament. Since the region r < r{t) charge-separated, 
the electric field due to the upstream particles at r = r(f) is 
Er ^ 2TTqnr and the magnetic field is Bg ~ lir/S^qrir. Let /i'(x) 
and /(x) denote, respectively, the linear particle number and 
the particle current density of the counterstreaming particles. 



— / n'2TTrdr, 



J 



cP'^nlirrdr. 



(4) 



Then the electromagnetic field of the counterstream is given 
by E'r ~ 2q'ii'/X and B'g ~ Iq'f /cX. 

Parameterizing r(f ) in terms of x = x(f ) and substituting the 
values of the fields in equation (|3} and dividing by mc^ we 
obtain 



dx 



dx I 



-27r(l — (3x)nrq — 2 



c r 



(5) 

The initial condition is r(-oo) = and we are interested in the 
evolution of r(x) up to the point where r = A. 

Particle number conservation for upstream particles im- 
plies that /J-ooH-oo = A"- Also note that qcl = -e^. Let 
= f lc[i! < denote the average velocity of the counter- 
streaming particles. Then the upstream particles experience 



^Ux ( 



7 



dr 



7- 



dx 



1 

2P 



(1- 



A). (1- 

r-\ 



(6) 



where we have defined the plasma skin depth of the upstream 
particles 5 = c/(47rn_ooe^/7-oo'M)'''^ and have divided the 
equation by 7-00 . 

Note that (j)(x, r) vanishes as r approaches A. Therefore it is 
a reasonable approximation to set the potential at the bound- 
ary of the charge-separated region to zero, which implies that 
7 ^ 7_oo is a constant and the upstream remains ultrarela- 
tivistic, X ~ P-ao ^1- "L&t X = x/5, R = f/6, L = X/5, 
and A^' = /i'/7r(S^n_oo; note that the latter quantity is equal to 
7_oo/4 times the so-called Budker's parameter. Then equation 
(jSJl can be rewritten as 

d^R ..... 



dx^ 



2L 



with the initial value R = Q and final value R ^ L. 

Equation Q describes the mechanics of charge separation 
among the upstream particles. The charge separation, and 
thus the self -pinching of the current filament, is possible only 
if the right hand side of the equation is positive, or 



N' > 



(l-A) 



RL. 



(8) 



Evidently, without the presence of a small amount of coun- 
terstreaming particles, the current does not pinch, and the fil- 
ament does not form. Equation Q can be solved once the 
density of counterstreaming particles, A^', is self-consistently 
specified. 

We use a toy model to integrate the trajectory of an up- 
stream particle. The electrostatic potential and the magnetic 
field are taken to be 



(j){x,r) = C^-, 
C7- 



./3 



Bz(x,r) = - 



2X 



(l + e.v/A)2 

(1+W^)2 



H(-x)+H(x) 



H(-x) + H(x) 



sm 



(9) 



with A = AX, 13 = 0.9, 7_oo = 20, and C = 0.85 is a factor tak- 
ing into account that final charge separation is not complete. 
Representative trajectories are shown in Figure |2] Note that 
the particle with initial axial distance r(-(X)) = 1.98A crosses 
r = A near x = 0, where charge separation must be maximum. 

6. THE COUNTERSTREAM 

The phase space distribution of the counterstream particles 
depends on the dynamics of the downstream x > xq where the 
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Fig. 2. — Trajectory of three representative upstream particles in the elec- 
tromagnetic field of the toy model (see text). Initial radial displacements of 
the particles are r(—oo) = (1, 1.8, 1.98)A. The grayscale shows the magnetic 
field, which vanishes atx = -oo and r = 0. The radial coordinate was extended 
to negative values to cover the meridional plane. 



u 




Fig. 3. — Trajectory of a counterstream particle in the electric field of 
the toy model (see text). The grayscale shows the electrostatic potential, 
which vanishes at x = -oo and r = it A. The particle has initial Lorentz fac- 
tor 7(-l-oo) ~ 15 and initial perpendicular velocity when crossing r = of 
/3, (— oo) = 0.3. The radial coordinate was extended to negative values to cover 
the meridional plane. 



electrostatic approximation breaks down. Therefore, the dis- 
tribution cannot be derived self-consistently within the scope 
of the present work. Instead, one can attempt to guess the 
density of the counterstream particles as a function of posi- 
tion. The most natural choice is that in which the density /i' 
is a function of the electrostatic potential only, because 0(jc) 
determines how many counterstream particles are energetic 
enough to reach x. There is not a unique distribution giving 
rise to a specific potential. The simplest density profile moti- 
vated by Maxwellian-like distributions is oc e"* "^/t"™'' , 
where 70 is a constant. 

Since 0(-oo) = 0, an exponential density profile implies that 
a finite number of unbound counterstreaming particles exist 
in the upstream. We here select only the particles that are 
confined by the electrostatic potential. We will return to the 
particles belonging to the high-energy tail of the Maxwellian, 
which are not confined, in § |8] Some particles can also be 
magnetically confined to a part of the filament even if they are 
not electrostatically confined. The particles that are not con- 
fined are not charge separated into filaments. Linear density 
of counterstreaming particles is given by 



^'(ji:) = ^'(ji:o)- 



-q'4>(x)/'romc^ . 



1 



-q'cl>(xo)homc- _ \ ' 



where (pix) is the electrostatic potential at {x,r ■■ 
depends on the charge density via 



(10) 
0), which 



q <p q 

mc^ mc 

= 2zr 
2 



^{TTr\2\n(X/r)+l]nq + 2n'q'} 



+N' 



(11) 



The density of the counterstreaming particles and the electro- 
static potential are transcendental solutions of equations J10> 
andlinil- 

The situation simplifies greatly if A^' <C and we can ig- 
nore the contribution of the counterstreaming particles to the 
potential. Of course, we must retain the contribution of the 
counterstreaming particles to the magnetic field in equation 
Q. To simplify this equation, note that deep in the up- 
stream R<^ L, while l3_, 1 and (1 - 1 Since 
1 < (1 -/3J.) < 2, to an order of magnitude we can replace this 
quantity with 2. Thus, after substituting equations ilQ\ and 
(^3 in equation Q and keeping the lowest-order terms in R 



we obtain 
(fR 



R 



■ + N'(xo) 



L \R 



-7-00 /."/470 



. (12) 



472^ ■ ' ' 270 

Current pinc hing is possible only when the right hand side 
of equation (I12t is positive — a condition that is manifestly 
not satisfied for very small R, roughly when R ^ 1/7^0^. 
Roots of this apparent paradox lie in the approximate treat- 
ment of the counterstreaming particles. We believe that the 
paradox would be eliminated by correctly treating the full 
three-dimensional orbits of the particles. 

In the regime in which the second term on the right hand 
side of equation M2\ is much larger than the first term, we 
estimate the width of the charge separation layer A by replac- 
ing the second derivative with a ratio, cf-R j dX^ ^ R/ A"^. If 70 
is interpreted as the temperature of the downstream medium, 
conservation laws imply that 7-00/70 ^ 3\/2- The width of 
the charge separation layer in units of the plasma skin depth 
then reads 



A. 



0.7 



R 




z'" 









-1-1/2 



0.53 L- 



(13) 



Note that the number of counterstream particles that are 
not confined by the electrostatic potential is given by A^^^.^^ ~ 

N'(xo)e~'^°°^ l'''^'^ . Therefore the square of the exponential fac- 
tor in equation ( I13> must be much greater than unity to pre- 
vent mass exodus of the counterstream particles into the up- 
stream. We will see in § however, that L itself cannot be 
much larger than unity or the Alfven limit would be violated. 
In §|8]we will find that stability requires N'(x^)) <C 1 in shocks 
with very large Lorentz factors. Therefore, the transition layer 
is many plasma skin depths wide. This underscores the chal- 
lenges that must be overcome to reproduce the steady-state 
transition layers in PIC simulations, which must resolve very 
many skin depths in the parallel direction. 

In Figure|3] we plot the trajectory of a counterstream parti- 
cle in the toy-model electromagnetic field of equation (|9}. The 
particle has initial Lorentz factor -fi+oo) w 15, just sufficient 
to keep it electrostatically confined. The initial perpendicular 
velocity of the particle when it crosses r = is Pri-oo) = 0.3. 
Note that the particle reverses at x ~ -15A. 

7. FILAMENT SIZE 

We have so far assumed that the radius A of the first genera- 
tion of filaments is constant throughout the charge separation 
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layer. In the standard picture motivated by cold shell colli- 
sions, filaments first form with radii on scales of the fastest 
growing mode of the transverse Weibel instability. The fila- 
ments immediately proceed to merge with each other in hier- 
archical fashion. 

In the collision of two cold shells of comparable density, the 
dispersio n rel ation for the filamentation mode derived from 

reads 



equation (IC5> by taking the limit — 

u'^-cViiO^ + Lul/j) 



■0, 



(14) 



where ujp = (Aire^n/me)^!^ is the plasma frequency written in 
terms of the total particle density in the shock frame n. This 
is an ultrarelativistic analogue of the nonrelativistic relation 
derived by Fried ( 19591). 

Equation MA\ implies that there is a growing mode asso- 
ciated with every wavelength. The growth rate decays on 
large scales lm{uS) oc A:'/^, and saturates at the constant value 
Im(ijj) ~ LOpl^^l'^ in the limit A: ^ oo. Although, formally, the 
infinitesimally small scales grow the fastest, the growth rate is 
nearly uniform on all scales larger than the plasma skin depth, 
i.e., for k > 2tt/6 where S = j^^^c/cup. Therefore, filaments 
with radius A > (5/4 form directly in the linear development 
of the transverse Weibel instability in cold shell collisions. 

A powerful argument, referring to the steady-state structure 
of the transition layer, argues for filamentation at the scale 
of the plasma skin depth. We showed in § |6l that the non- 
neutrality of the current filaments gives rise to an electrostatic 
potential that keeps most of the thermalized particles from es- 
caping into the upstream domain. For the confinement to be 
possible, the depth of the potential well \(j>\ ~ TrA^xne must 
be larger than the temperature of the downstream medium 
~ 7o'«c^, where x ^ 1 is the fraction of particles confined 
to a filament when charge separation is maximum and 70 ^ 
7_oo/3V2- Therefore, A > 2x~^^^S, which sets lower limit on 
the filament radius to be about twice the plasma skin depth. 

On the other hand, the current inside the filament |/| ~ 
A^X«^c cannot exceed the Alfven limit 7^ ^ A-fmc^ /e, where 
A ^ 1 is a numerical form factor. (In principle, the current can 
exceed the limit if the perpendicular velocities in the filament 
greatly exceed the parallel velocities, which is not the case 
here.) From this one can conclude that A < 

2^1/2^-1/2^^ Evi- 

dently, the limit is saturated when charge separation is almost 
complete in a filament with radius just slightly larger than the 
plasma skin depth. The above two constraints pen the radius 
of the filament to be just slightly larger than the plasma skin 
depth, A ^ 26, where we have taken into account that practi- 
cally all particles become confined within filaments, x ~ 1, at 
the end of the charge separation phase. 

Similarly, for the magnetic field |B| ^ 2\I\/\c to be dy- 
namically important, the ratio of the magnetic to mechan- 
ical and thermal energy density = B- /^Trjmc^n should 
be of the order of unity. This ratio can be expressed as 
eB 3|///aP(A/(5)-2. sjjj^g < follows that A > 
S would imply ^ 1, rendering the magnetic field dy- 
namically unimportant. Indeed, ^ 0.1 - 1 has been ob- 
serve d in PIC simulations of shell collisions in pair plasmas 
(e.g..lKa7.imura et al.l1 998tln;^ov>2nnk.Silva et al■l20n:l^ 



iMedvedev etalJ2005l) 



8. THE PRECURSOR 

If the counterstreaming particles at xq are thermal with 
a Maxwellian tail extending to energies j'mc^ > q'(l)(x()), a 
steady stream of particles escape into the upstream. These 



particles run ahead of the shock wave as a "precursor." In 
principle, the precursor can initiate transverse Weibel insta- 
bility and thus transfer energy ahead of the shock. This would 
invalidate the premise that the upstream is cold and field-free 
as it enters the charge separation layer Fortunately, finite per- 
pendicular motions of the counterstream particles in the up- 
stream frame can suppress the instability if the density of the 
escaping particles is not too large. We here discuss conditions 
for the stability of the precursor 

The approximate, waterbag-type dispersion relation for 
transverse electromagnetic modes, assuming that the up- 
stream is cold and the precursor is warm, is derived in Ap- 
pendix |H1 To check whether the upstream is stable in the 
presence of a low density precursor, we consider the case 
A^fi-ee ^ 7-00 where A^^ee = n'/n-oo is the ratio of the density 
of the escaping precursor to the density of the upstream in 
the frame of the shock. With this assumption, the upstream 
is nonrelativistic in the special frame in which the perpendic- 
ular electric field vanishes (see App endi x 151. The unstable 
wave numbers must satisfy equation ( IC3t and thus instability 
takes place when (an equivalent relation was also derived in 
iSilva et alJ20()l 



-00/7-c 



2 



(15) 



where f3'j_ is the typical perpendicular velocity of the pre- 
cursor particles in the frame of the current defined in § |3 
which is related to the same quantity in the shock frame via 
= l3'Jj(l + < 1/7-00- Given that majority of the 
precursor particles have 7' ~ 7-00^ the upstream-precursor 
system is stable as long as N^^.^^ < 7!^ and is unstable for 
^free ^ 7-00 ■ Therefore to check whether the precursor initi- 
ates transverse Weibel instability ahead of the shock, we must 
compare A^fj.gg with 7!^. 

The requirement that the precursor be stable to the develop- 
ment of the transverse Weibel instability can be seen as set- 
ting a condition on the number of particles that are contained 
in the unconfined, high-energy Maxwellian tail. Within the 
model presented in §|6l the density of these particles is 



-L-7_oo/47o 



(16) 



In an ultrarelativistic shock, energy conservation dictates that 
70 ^ 7_oo/3 V2- We have argued in §0that L^2. Then from 
equation (O, we have that Nf^^jN'ixo) - e"^^^ « 0.014. The 
precursor is stable if Nf^.^^ < ^l^, or ^I^N'{xq) < 70. 

Thus if 7_oo > e^/^ ~ 8, the precursor can be stable only 
if A^' < 1 , which means that the density of the counterstream 
particles at xq is less than the density of the upstream. How- 
ever the density of the counterstreaming particles must be- 
come comparable to that of the upstream at some point in the 
transition layer, since one-third of the downstream particles 
have negative velocities in the frame of the shock. There- 
fore in shocks with large Lorentz factors, compression and 
isotropization must take place downstream of the charge sep- 
aration layer 

We speculate that if the upstream medium is weakly magne- 
tized, reflection of precursor particles back toward the shock 
might play a role in the injection of particles into the diffusive 
accelerator mechanism of the shock front. 

9. DISCUSSION 
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In the hydrodynamic limit, the flow crossing an adia- 
batic relativistic shock wave is compressed by the factor 
n+oo/n-oo ~ 3. We saw that the structure of the transition 
layer is sensitive to the strength of the electrostatic poten- 
tial and the width of the current filaments. If the potential is 
strong enough to decelerate the incoming upstream to a non- 
relativistic mean bulk velocity in the frame of the shock (note 
that each particle is still relativistic), compression can take 
place in the electrostatic charge separation layer, although, of 
course, thermalization still requires a dynamical electromag- 
netic field. If the flow remains relativistic through the charge 
separation layer, the shock compression must take place after- 
wards, and the charge separation layer is followed by a com- 
pression layer We elaborate briefly on the latter possibility. 

In the charge separation layer, current filaments do not in- 
teract because the perpendicular displacement of the upstream 
particles is of the order of the filament size. Also, there is a 
near-perfect balance of the magnetic field (which attracts like 
currents) and the electric field (which repels them). Therefore, 
the charge separation layer is approximately two-dimensional. 
This two-dimensional picture breaks down when the upstream 
particles have crossed the filaments once and charge sepa- 
ration has reached maximum. At that point, the filaments 
start to interact and the picture becomes three-dimensional. 
These three-dimensional processes also govern the decay of 
the magnetic field. 

Magnetic interactions cause rearrangement among the fil- 
aments (unfolding, mergers, etc.). Even noninteracting fila- 
ments do not remain parallel to the direction of propagation 
of the shock. Based on an analogy with magnetohydrodynam- 
ical instabilities in magnetically-confined columns, initially 
straight coUisionless current filaments should be unstable to 
a kink or pinch-like mode. The filaments bend, buckle, and 
fold. Time-varying electric fields associated with the rear- 
rangement scatter particles inelastically to large pitch angles, 
thereby aiding diffusion of the momentum, rendering the flow 
isotropic and eventually thermalized in the frame of the down- 
stream plasma. Isotropization in the downstream frame im- 
plies compression. 

The results thus derived apply to e~e^ shocks. Does the 
same picture describe electron-ion shocks? Simulations of 
colliding electron-ion shells ( Fred eriksen et al. 2004) reveal 
that the electrons thermalize before the ions do. If the transi- 
tion layer contains ion filaments, as the simulations suggest to 
be the case, the electrons are scattered by the electromagnetic 
field of the ion filaments. Charge separation among the ions 
takes place in the background of hot electrons. The mecha- 
nism of charge separation presented in §|5]and §|6lcould ap- 
ply to the ustream ions. They pinch into filaments as a small 
amount of hot electrons are trapped within each upstream ion 
filament and cancel some of the ionic self-electric field. The 
upstream ion filaments are much like the e~e^ filaments be- 



cause the downstream electrons moving toward the upstream 
have comparable energies, and thus comparable inertia, as the 
downstream ions. At this point we cannot explain the force 
that keeps the downstream ions from escaping into the up- 
stream, although we expect that an electrostatic field associ- 
ated with the upstream electrons is involved. 

10. CONCLUSIONS 

We here summarize the main conclusions of the paper: 

1 . The first generation of magnetically-confined current fil- 
aments that arise in the saturated state of the transverse Weibel 
instability in coUisionless shocks are nonneutral and nonuni- 
form. 

2. Two interpenetrating components in the shock transition 
layer are the upstream particles and the downstream particles 
moving toward the upstream in the frame of the shock ("the 
counterstream"). Both components are charge-separated. 

3. Space charge in the filaments gives rise to an electro- 
static potential; we call the leading, electrostatic edge of the 
transition layer "the charge separation layer." 

4. In the charge separation layer the charge density and the 
magnetic energy density increase in the parallel direction. 

5. The electrostatic potential confines most of the counter- 
stream particles to the charge separation layer, thereby pre- 
venting their escape into the upstream. 

6. The radii of first-generation filaments must be just 
slightly larger than the plasma skin depth. 

7. A small fraction of downstream particles, "the precur- 
sor," escape into the upstream. These particles may be impor- 
tant for diffusive acceleration in the shock. 

8. The precursor moving through the upstream must be sta- 
ble to the development of the transverse Weibel instability, 
which sets a condition on the density of counterstreaming par- 
ticles in the charge separation layer 

9. Shock compression is completed downstream of the 
charge separation layer, rather than in the layer itself. 

10. Downstream of the charge separation layer filaments 
merge on a dynamical time scale and the electromagnetic field 
is dynamic, facilitating isotropization and thermalization of 
the plasma. 
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APPENDIX 

DIELECTRIC PERMITIVITY TENSOR 

We work in the orthonormal basis (e||,ej^i,ej^2) where e|| is parallel to the direction of propagation of the shock. Consider 
linear electromagnetic perturbations of the field-free plasma of the form 

E(x, t) = E(k, u) exp[;(k • x - cot)] , 

B(x, f ) = B(k, uj) exp[i{k ■ x - ujt)] , (Al) 

subject to the condition k x E(k,cj) = (w/c)B(k,cj). To study filamentation in the transition layer of the shock, we consider 
transverse modes, k^^ = 0. This is a Lorentz-invariant choice as long as u; is purely imaginary. Components of the dielectric 
permitivity tensor e can be evaluated using standard methods (e.g., iKrall & Trivelpiecall973h . The perturbed electric field 
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satisfies the linearized Maxwell- Vlasov equations D • E = where E = (E^^,E±i,E±2) and D,y = eij + (c^/uj^){kikj-k^6ij). Note 
that £ and E refer to the permitivity tensor and the electric field in Fourier space; explicit dependence on the wave vector is 
omitted here. We use the index a to label distinct species of particles with phase space distribution functions fa(p), particle 
rest masses OTq, concentrations in the observer frame iia, and plasma frequencies uipa = (AimaC^ /niayi^ . Components of the 
permitivity tensor then equal 

,2 



-5, 



where the velocities are related to the momenta via p = jma\. 

DISPERSION RELATION 
We choose to work with the waterbag distribution function for each plasma component 

1 

TTpj_ 

where is the bulk momentum of the plasma, p±a is the maximum perpendicular momentum, 6{x) is the Dirac (5-function 



fJp) = ^:^S(p\\-p\\Mplo.-pil (Bi) 



and H(x) is the Heaviside step function. The distribution function in equation iBH assumes that the plasma is cold in the parallel 
direction. Note that t he dist ribution is isotropic around the a xis of propagati on of the shoc k. The w aterbag-type distribution func- 
tions in Medvedev & Loebllfl99 9). Silva et al. (2002), Fons eca et alJ I l20()3h . and Wiers ma & Achte rberg (2004) are not isotropic 
in this way. The differences stemming from the particular choice of distribution function, however, do not seem to affect the 
instability growth rates substantially. 

Current filaments parallel to direction of propagation of the shock form from electromagnetic modes with purely perpendicular 
magnetic field, B|| =0, and purely parallel electric field, Ej^ = 0. The latter condition is violated in a general reference frame. 
There exists a unique Lorentz frame in which the parallel and perpendicular electric fields decouple and the linearized Maxwell- 
Vlasov equation is block-diagonal, that is, the cross term vanishes = 0. In this frame, the dispersion relation for the linear 
development of the filamentation instability is simply = 0. Substituting equation dBlt into the first and second of equations 
( lA2t yields 



2(l+7a7||a-7p 2 

■+P\\a 



£)|lll=l ^ — 

Q 



where, omitting subscript a, we used the definitions = p\\/^mc, f3± = p^/^mc, (3"^ = /3| + /3^, 7 = (1-/3^) 'Z^, and 7|| = 
(l+7'/3|)'/'. 

Consider large phase velocity perturbations Lojk^ cf3± in cold plasmas 7/3 <C 1. Note that 
and 




ck\ 3 /ck 
oj 1 4 V w 



4 

Pl + 0[(ckp^/u)^] 



2lA 



c^k 



- + O[{ckp^/uo)\ (B4) 



Thus, to the lowest order in the condition j^,- = 0, required for Ei_ = and B|| = 0, reduces to 

E— /5|t"«E— /3|l" = 0, (B5) 

^ la ^ ma 

a a 

where we have defined rest-frame particle concentrations = na/la- The frame in which the perpendicular electric field 
vanishes is not the rest frame of the plasma. 

STABILITY 

Consider a plasma containing a dense, infinitely cold component, and a cou nterstreaming low-density component with finite 
perpendicular temperature. Following the notation of lAkhiezer et a we add a prime to the quantities pertaining to 
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the low-density component and dispense with the subscript a. Assume that nf-f ^ n'/7'- Then the low-density component 



is ultrarelativistic, 7' ^ 1 and /Jf, ^ -1. Equation (IB5> then implies that uj^jS^^/'-f ^ uj'J'/j'. Defining the ratio of the rest- 



frame concentrations a = n[JnQ <C 1, we have that /3|| ~ a <C 1 and 7 ~ 1. In the same context. FSilva et all J2002h define the 

concentration ratio in the frame in which the denser component is at rest. 

In the frame with a vanishing perpendicular electric field, the dispersion relation simplified using JB4> and reads 



or, using ui' /"f' ~ aWp and keeping the lowest order terms 



l + ( - 



- , + T =0. (CI) 



p ' ' p 



l-( — ^ r=0- (C2) 

This equation is stable for all values of the wave vector when 4a < 3/3^^. For 4a > 3/3^^, the condition for instabihty is 

The growth rate and the wave number of the fastest growing mode are 

^L. « I \l t/?1 I , « ^ I a/ - 1 I • (C4) 




Next, consider the case in which both components of the plasma are relativistic, have similar densities, and one component is 



infinitely cold while the other has finite perpendicular temperature. Here /3|| ^ (3',, ^ I and thus ujI/"/ ^ ti^' /Y- Relevant terms 



in the dispersion relation then read 



The condition for instability is then 



9 ^ 16 ujI 

k^<- (C6) 

The growth rate and the wave number of the fastest growing mode are 

(C7) 

Note that here ujp and 7 are the plasma frequency and Lorentz factor of the infinitely cold component. We made no assumption 
about which component is denser, but have assumed that the densities are similar enough for both components to be relativistic 
in the frame in which the perpendicular electric field vanishes. 
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